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Abstract 

We investigate the amount of fine tuning of the electroweak scale in the presence of new 
physics beyond the MSSM, parametrized by higher dimensional operators. We show that 
these significantly reduce the MSSM fine tuning to A < 10 for a Higgs mass between the 
LEPII bound and 130 GeV, and a corresponding scale M* of new physics as high as 30 
to 65 times the Higgsino mass. If the fine-tuning criterion is indeed of physical relevance, 
the findings indicate the presence of new physics in the form of new states of mass of 
0(M*) that generated the effective operators in the first instance. At small tan/3 these 
states can be a gauge singlet or a SU(2) triplet. We derive analytical results for the EW 
scale fine-tuning for the MSSM with higher dimensional operators, including the quantum 
corrections which are also applicable to the pure MSSM case in the limit the coefficients 
of the higher dimension operators vanish. A general expression for the fine-tuning is also 
obtained for an arbitrary two-Higgs doublet potential. 
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1 Introduction. 



Low-energy super symmetry offers an elegant solution to the hierarchy problem. One conse- 
quence of this is that it introduces a spectrum of supersymmetric states in the visible sector 
with mass of the order of the electroweak scale. However, none of the superpartners of the 
Standard Model have been seen, although there is hope that LHC will soon remedy this. In 
trying to determine the physics beyond the Standard Model the fact that no superpartners 
have been observed is significant as it (re) introduces the need for some amount of fine tun- 
ing of the parameters of the minimal supersymmetric standard model (MSSM), to separate 
the electroweak and supersymmetry breaking scales (the "little hierarchy problem" ) . On the 
other hand circumstantial evidence for supersymmetry such as the successful unification of 
couplings [U O O Hj or radiative electroweak breaking [3] is consistent with and in fact requires 
the existence of such light superpartners. 

The basic issue raised by fine-tuning is the sensitivity of the electroweak scale (more 
precisely the mass of Z) to small variations of the input parameters of the MSSM, consistent 
with the measured Z mass, mz, and the current bounds on the lightest Higgs mass, mn- The 
need to keep fine tuning small indicates a light Higgs in some tension with the LEPII bound 
[Sj m h > 114.4 GeV. Consistency of this bound with the MSSM tree level bound rrih < mz 
can only be achieved at the quantum level, by a large top quark/squarks loop correction to 
mh- To maximise this, the top squarks must be quite massive or highly mixed implying that 
the MSSM model is fine-tuned. 

Various definitions of fine tuning have been proposed. The most popular one [7] is based on 
the logarithmic derivatives of the observables with respect to the set of parameters considered. 
It has been widely used in quantifying the fine tuning in the MSSM [8} 19] flU l fTTl \12\ fT3l fli] 
and we shall use it in this analysis. The new feature of our analysis starts with the premise 
that if low-energy supersymmetry is indeed the solution to the hierarchy problem, a significant 
amount of fine tuning of the electroweak scale in the MSSM may in fact suggest that there are 
additional new degrees of freedom in the theory beyond those of the MSSM. There are many 
models that consider additional degrees of freedom beyond the MSSM in order to reduce the 
amount of fine tuning. The NMSSM is just such an example [15] which has an extra chiral 
singlet. One can consider other MSSM extensions with more chiral fields, additional gauge 
interactions, etc. Each of these brings different solutions to the little hierarchy problem and it 
is difficult to assess which of these is the most compelling. In this paper we perform a model 
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independent analysis of the nature of this new physics based on a general parametrisation 
of physics beyond the MSSM. In particular we extend the MSSM by the addition of higher 
dimensional operators [HI [TTJ, [18j [19] that encode the effect of all possible new physics at 
scales below the appearance of the new degrees of freedom. Having identified the most relevant 
operators one can later address the question of what new physics generated these operators 
in the first instance. The advantage of the effective approach is that it provides an organising 
principle according to which one usually restricts the analysis to operators of a given (leading) 
order in the scale of new physics M*, with higher order operators suppressed by higher powers 
of M*. The analysis we consider includes dimension d = 5 and d = 6 operators beyond MSSM 
[20\ [2T| 123] . For the case of d = 5 operators we determine the amount of fine tuning as a 
function of the mass of the lightest Higgs corrected by the quantum contributions using both 
analytical and numerical techniques. One-loop renormalisation group corrections in the Higgs 
potential are also included. 

The expectation that higher dimensional operators can reduce the amount of fine tuning is 
broadly based on two arguments. Firstly these operators may directly increase the tree level 
value of to/j |20l l21j. Consequently the tree level upper bound on may be relaxed, and the 
quantum effects needed to satisfy LEPII bound may be smaller, corresponding to reduced fine 
tuning. Secondly, the higher dimensional operators may generate additional contributions to 
the quartic Higgs couplings of the MSSM, again serving to increase the Higgs mass. This 
effect can be quite significant because, in the MSSM, the quartic coupling, (g% + gf)/8, where 
52 and gi are the SU(2) x U(l) gauge couplings, is anomalously small; indeed its smallness 
is a major source of the little hierarchy problem. For the case of just the d = 5 operators 
a numerical study shows that these effects can reduce the amount of fine tuning, A, of the 
electroweak (EW) scale relative to the MSSM case, to less than 10 for a Higgs mass in the 
range 114.4 GeV < < 130 GeV even for a scale of new physics as high as (30 to 65) times 
the higgsino mass, and possibly above the LHC reach. We also give in Appendix an analytical 
formula for the EW fine-tuning in a general two-Higgs doublet model, which can be easily 
applied to specific models. 

The plan of the paper is as follows. Section [2] lists the d = 5 and d = 6 operators that 
are consistent with the MSSM symmetries and that can affect fine tuning. In Section [3] we 
evaluate analytically and numerically the fine tuning in the MSSM extended by the d = 5 
operators. The conclusions are given in Section HI 
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2 Higher dimensional operators beyond MSSM Higgs sector 



In this section we list the effective operators of dimension d = 5, 6 that can be present in the 
Higgs sector consistent with the symmetries of the MSSM. These operators parametrise new 
physics beyond the MSSM and affect the Higgs scalar potential. Therefore they also affect 
the amount of fine tuning of the EW scale, as discussed in detail in the next section. The 
(R— parity conserving) d = 5 operators in the MSSM Higgs sector are: 

C 2 = J d A 9 { A(S,S^)D a B(S,S^)H 2 e- Vl D a C(S, 5*) e Vl H x +fc.c.} (2) 

where S is the spurion field, S = 99mo, A(S,S^), B(S,S^), C(S,S^) are polynomials in 
S, and niQ is the susy breaking scale in the visible sector (in gravity mediation mo = 
{Fh) /Mpi anc k where (F^) is the auxiliary field vacuum expectation value (vev) in the hidden 
sector responsible for supersymmetry breaking). As we discuss in Section [331 the first operator 
can be generated, for example, by integrating out massive gauge singlets or SU (2) triplets, 
while the second is easily generated by integrating out a pair of massive Higgs doublets [21) . 
all of mass of order M* . 

In [21[ [22] it was shown that by using general field redefinitions one can remove C 2 from 
the action. The effect of this is an overall renormalisation of the soft terms and of the /j, 
term. Since the fine tuning measure includes the fine tuning with respect to each of these soft 
operators separately adding C 2 cannot reduce the overall fine tuning. For this reason we will 
only include C\ in our discussion of fine tuning with d = 5 operators. 

There are also d = 6 operators that can be present in addition to the MSSM Higgs sector. 
These are suppressed relative to the d = 5 operators by the factor 1/M*. However they may 
give contributions to the Higgs potential enhanced by tan/3 relative to the d = 5 so cannot 
be ignored at very large tan/3. The list of d = 6 operators is (see also [161 02] ) : 

O l = -^J^eZi^^iHje^Hi) 2 , i = l,2. 

° 3 = M I dH Zz{S) s]) {H ' eVl Hl) {H " eV2 H2} ' (3) 

(These can be generated by integrating a massive £7(1) gauge boson or a ££7(2) triplet). 
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° 5 = wj di ° Z5(5, 5t) (jffit eV/l Fl) ( ^ 2 ^ + /i - c - ) 

O7 = ^ J d 2 6 Z 7 (S,0)W a W a (H 2 H 1 ) + h.c, 

O s = JflJ^ [M0,&) {H^Htf + h.c^ (4) 

where W a is the super symmetric field strength of a vector superfield of the SM gauge group. 
04 can for example be generated by integrating a gauge singlet. 

Oio = ^2 y ^ ^io(5, 5+) F 2 f V 2 e^ 2 V 2 ff 2 

O12 = J d 4 6 Z l2 (S,tf) Hle V2 V a W a H 2 (5) 

where V Q acts on everything to the right and V a -ffj = e~ Vi D a e Vi Hi. i=l,2. In addition to 
the spurion dependence in the wavefunctions Zi(S, Sr) , extra (S, S 1 ') dependence (not shown) 
can be present under each derivative V a in eq.©, in order to ensure the most general super- 
symmetry breaking contribution associated to these operators. One may use the equations 



Note 



of motion to replace the operators involving extra derivatives by non-derivative ones 
that when computing the fine tuning measure eliminating a particular operator will lead to 
correlations between the remaining operators that, strictly, should be taken into account. 

Given the large number of d = 6 operators, determination of the fine tuning with respect 
to their coefficients is difficult. For this reason we restrict the following discussion to the d = 5 
operators. In Section 13.31 we comment on the new contributions that may come from d = 6 
operators and discuss the limit on our analysis that follows from keeping only d = 5 operators. 



Setting higher derivative operators onshell is a subtle issue in this case. One can also use general spurion- 
dependent field redefinitions to "gauge away" (some of) these operators, using the method of |21l 122] , 
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3 Fine-tuning in MSSM with d=5 operators (MSSM 5 ). 



3.1 The scalar potential 



In this section we evaluate the EW scale fine-tuning in the MSSM extended by C\ in eq.([T]). 
Including it together with the MSSM, the full Higgs Lagrangian is then given by 



rf 4 



l-c hl S S*) H\ e Vl H x + (l- c h2 S S ] ) H ] 2 e V2 K 



The corresponding scalar potential is given by 



/i (1 + B S) ffiff 2 + -^-(l + c S) (H! H 2 f 



+ h.c. 



(6) 



2 2 

V = ml\h 1 \ 2 + ml\h 2 \ 2 -(mlh l h 2 + h.c) + 9 —(\h l \ 2 -\h 2 \ 2 f + 9 —5\h 2 \ A (7) 

o o 

+ (|/ii| 2 + | h 2 | 2 ) (Ci hi h 2 + h.c.) + l - ( C 2 (hi h 2 ) 2 + h.c) 

2 — „2 _i_ „2 



where g 2 = gf + g$ 



and 



m?(i) 
m|(t) 
m§(i) 



Ci = 2 Hq/M* , C2 = -2 c m /M* (8) 
m + Mo °f(*) + m 12 ci(*) 

Mo^iC*) + m 2 l2 <Ji(t) +A t m mi2<7 5 (i) + m 2 o- 7 (i) - m 2 ,^ 2 <x 6 (i) 

Mo mi2 cr 2 (t) + -Bo m-o Mo o"8 (*) + Mo mo -4* 03 (i) (9) 



The coefficients <7j depend on t = lnM^/Q 2 with functional dependence given in [71 
[271 128] . The (high scale) boundary values (t = 0) are normally chosen to be o"i,2,..,6 = 0, 
°"7,8 = 1 (i-e. Ch\ 2 = 1 )• F° r Q 2 = m z (* = the values of these coefficients are given in 
Appendix [X] in terms of the top Yukawa coupling. To simplify notation we will not display 
the argument t z in what follows. 

The quartic term <5|/i 2 | 4 is generated radiatively |131 124] , Including leading log two-loop 
effects one has 



3/4 
g 2 it 2 



, Mr X t 

In — ^ + — + 
m t 4 32tt 2 



1 / 3/ i 2 -16 5 2 )(x t + 21n^)ln^ 
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= 2 (A t m -/icot/3) 2 / (A t m - fi cot /3) 2 \ 
1 M| V 12 M- 2 / ' 1 j 

with M~ = mj 2 , and 33 is the strong coupling. 
The minimum conditions for V can be written as: 

9 m 2 / , dm 2 9 <9A \ . , 

with the notation v 2 = v \ + v 2 , tan /3 = V2/V1, m 2 z = g 2 v 2 /4 and where: 

m 2 = rh 2 cos 2 /3 + m 2 ; srn2 P ~ m i s i n 2/3 

2 /• 

A = (cos 2 2/3 + 5 sin 4 /3) + Ci sin 2/3 + sin 2 2/3. (12) 
8 4 

Note that in deriving these expressions we have discarded non-leading log corrections except 
those to the quartic Higgs coupling where the tree level term is anomalously small. 

3.2 Analytical results for fine-tuning 

The fine tuning of the EW scale with respect to a set of parameters p introduced in [7] is 

d In v 2 

A = maxAbs[AJ , A„ = (13) 

With m 2 = m 2 (p,f3), X = X(p, j3) we can find dfi/dp from the second of eqs. fjlip (more 
precisely this determines the parameter dependence of /3 m i n ): 

dp If ^dXdm 2 2X^ 2m2 ^ m2 9X 2 \ 
dp z \ dp 9/3 dj3dp dp df3 d/3dp J 



where 



d 2 m 2 . 2 <9 2 A\ v 2 fdX\ 



Using this one finds 

f d 2 m 2 9 d 2 X\ f dX 1 dm 2 \ dm 2 d 2 X <9A<9 2 m 2 l . . 

+ V 2 — )(— + — — ] +— — T^TT- ■ (14) 



A p = - P - 

z 



d(3 2 df3 2 J \ dp v 2 dp J d(3 djidp df3 djidp 
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3.2.1 A general two-Higgs model 



Using eq. (|14p we derived a general analytical result for the fine-tuning of the EW scale in a 
general two-Higgs doublet model allowing for the most general renormalisable Higgs potential, 
see Appendix[Bl eq. (|B-6p . The results are presented in terms of derivatives of the soft masses 
and couplings of the scalar potential. 



3.2.2 The MSSM with dimension-five operators (MSSM 5 ) 

Applied to the case of the MSSM with dimension-five operators the results presented in 
Appendix [B] give: 



Mo 



hi' 



v 2 cos 2/3 sin 2/3 (1(272 -5gfv 2 /8) 



+ 2 7 i [5g 2 /8- ((4 + 5) 5 2 /8-C 2 ) cos2/3 
74 - v 2 (2Ci sin 2/3 - C2 cos 4/3)]} 



+ 2 



2/igtrf + (G ^-71) sin 2/3 



(15) 



■ -rvr; \ ~ A2 sin 4/3 [471 cos 2/3 + (5 g 2 v 2 /8 - 2 72 ) sin 2/3 



+ 2v 2 



2 (71 - // mi 2 cr 2 ) cos 2/3 + 73 sin 2/3 



4Ci cos 2/3 +5g 2 cos /3 sin 3 /3 



+ (C 2 -5 2 /2) sin 4/3 



74 - (2Ci sin 2/3 - C2 cos 4/3) 



O O 000/ I 1 

2m - 73 sin p + C2 « sin /3 cos /3 — (71 — m i2 /iq 02) sin 2/3 ^ 



(16) 



A., 



mi2 r v z 
I "2 



2/^oO"2 cos 2/3 — (j4 4 (T5 mo + 2 777,12 (0-4 — <ti)) sin 2/3 



4d cos 2/3 + 5 g 2 cos /3sin 3 /3 + (C 2 -5 2 /2) sin 4/3 +2 2 77112 01 - /-to <72 sin 2/3 
+ (A t cr 5 m + 2 mi 2 (ff4-ffi)) sin 2 p 74 - t; 2 (2d sin2/3 - C2 cos4/3) } (17) 
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At 



v 2 D 

+ 2 Civ 2 sin 2/3 -74 



1 2 niQ sin /3 2fiQ 03 cos /3 + (2 erg mo — 05 77742) sin/3 — (2 ^ 2 cos 4/3 



+ mo u 



//0C3 cos 2/3 — (1/2)05 mi2 sin 2/3 + 06 mo sin 2/3 



x 



4Cicos2/3 + <5c/ 2 cos/3sin 3 /3+ (C2 -5 2 / 2 ) sin4 /3 } 



(18) 



A R2 = - 



Also: 



2B °™^ a * {v 2 [2 Cl + (C 2 - (4 + S)g 2 ,8) sin 3 2/3 

{5 g 2 v 2 /16 - 72) sin 4/3 - 4 7i sin 2 2/3} 



+ 



^i 2 



2D' 



C2 



8L>' 



2 m§ - 2 cos 4/3 (3ml + (4 + 5) (0 V/8) sin 2/3) 

+ (3 (ml - m 2 ) + 5g 2 v 2 /8) sin 4/3 

7 sin 2/3 [ - 2 cos 4/3 (4 mj + (4 + 6) (g 2 v 2 /8) sin 2/3) 



+ ( 4 (ml - m\) + 5g 2 v 2 /8) sin 4/3 



with the notation: 



(19) 



(20) 



(21) 



D = 2 
- 2 



— w 
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4Ci cos 2/3 + C2 sin 4/3 + a 2 (5 cos /3 sin 3 0-1/2 sin 4/3) 



(22) 



Cisin2/3 + C 2 /4sin 2 2/3 + c/ 2 /8(cos 2 2/3 + 5sin 4 /3) v 2 (2d sin 2/3 - C2 cos 4/3) -74]} 



El 
4 



L>' = — (cos 2 2/3 + <5 sin 4 /3) (2 {m\ - fh\) + 5g 2 v 2 /8) cos 2/3 



- (4 + 5) (s V/8) cos 4/3 + 4 m^ sin 2/3 



2 „• 



5V 
32 



[2- (4 + 5) sin 2 /3] 2 sin 2 2/3 



(23) 
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and 

71 = fJL (B m a 8 + mi 2 2 + A t m a 3 ) 

72 = (-1 + <?7 - A% C7 6 ) ml + A t cr 5 m mi 2 + m\ 2 (04 -01) +5 # 2 v 2 /16 

73 = 2 (1 - 07 + A 2 <7 6 ) m 2 , - 05 mi2 m 

74 = 2 72 cos 2/3 + 471 sin 2/3 - (4 + 5) (g 2 v 2 /8) cos 4/3 (24) 

The contributions A^2 are proportional to Q so, for small enough changes from the MSSM 
case, the fine-tuning introduced with respect to these new parameters is small and sub-leading 
relative to that for the other parameters. 

It is convenient to treat /3 as the free parameter rather than Bq . Using the second minimum 
condition of ([11]) (after replacing m 2 by Q), one finds 



B = — {/x o mi202 + Mo rn Q A t a 3 - \{m\ + ml) 



sin 2/3 



+ z o ^ z 5 ( Ci (1 + sin 2 2/3) + ^ sin 2/3 + J (g 2 /8) sin 2/3 sin 2 /3 ) [> (25) 
mf + m| V 2 



Note that 7^4 brings some extra £1,2 dependence through £?0) while 72,3 are Ci,2 independent. 
Ap> p = {/-ig, mQ, j4 2 , Sq, m 2 2 } contain some 0{C,f) terms, although the potential is only linear 
in £j. In the above expressions all coefficients 01,2,. ..,8 are evaluated at mz and their values are 
given in (jA-ip . They depend only on the top Yukawa coupling at mz- The only approximation 
in obtaining the above expressions for A p 's is that we did not include the effect of derivatives 
(with respect to parameter p) acting on 5 (the radiative correction to the quartic term). This 
is a legitimate approximation since this effect is numerically very small (for the MSSM alone 
it induces an error for fine-tuning A equal to or less than unity, while in the MSSM5 the error 
is even smaller (1%); for larger tan/3 this error is further reduced). 

The above results for the fine-tuning measure simplify in the limit of ignoring the RG 
effects on the masses i.e. 01,2,3,. ..,6 = 0; 07,3 = l.In this case 

71 = fi B m , 72 = 5g 2 v 2 /l6, 73 = 0, 

74 = 4^0 7710^0 sin 2/3 + (g 2 v 2 /8)[ 5 cos 2/3 -(4 + 5) cos 4/3] (26) 

Ignoring RG effects on the quartic couplings too, 5 = 0, then 74 = 4Z?o mo M) sin 2/3 — 
(g 2 v 2 /2) cos 4/3, 72 = 0. Finally, in the limit Ci,2 = of the fine tuning relations A p , one 
obtains analytical expressions for the EW scale fine tuning in the MSSM alone, with 71 ...4 as 
in (j24"|) . Since these may be useful for other studies, they are provided in Appendix lAl 
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3.3 The large tan/3 limit. 



The above formulae for the fine-tuning simplify considerably in the limit of large tan f3. Ig- 
noring terms suppressed by inverse powers of tan (3 one has 



A „ 2 



A., 



A 2 

m 12 



A, 



(1 + 8) ml 

-m 
(1 + 5) m| 

-mi2 
(1 + 8) m% 

A t m 



2a 7 m + A t (a 5 m 12 - 2A t cr 6 mo) + 

Ci v 2 fj, cr 2 



/'o 



Cl V 2 77112 0"2 



m 72 + (1 + 8/4) m| 



^4t a 5 m + 2 <r 4 m\i - 
2 A t cr 6 m - (7 5 77112 + 



A 



p, 2 



\ 2 



(1 + 8) m 2 

-Ci ("712 Q-2 + A t m a 3 ) fj, v 2 
(1 + 5) m% (72 +m%(l + 6/ A)) 



72 + (1 + (5/4) m| 

Cl Mo CT3 

72 + (1 + 5/4) m| 



0, i 



1,2 



(27) 



One may see that all the fine tuning measures are suppressed by the factor (1 + 5) _1 demon- 
strating why quantum corrections to the quartic Higgs coupling can significantly reduce the 
fine-tuning. A similar effect applies at small tan f3 as well. 

3.3.1 Dimension-six operators 

The operator analysis used here has a limited range of validity because it corresponds to 
integrating out new heavy degrees of freedom. If the mass of these degrees of freedom is not 
much above the energies being probed, the operator analysis breaks down and one must deal 
with the new degrees of freedom directly. The mass, M*, at which this happens corresponds 
to the point where high dimension operators are not suppressed relative to low dimension 
operators. A measure of this may be obtained by dimensional analysis in which the operator 
matrix elements are taken to be determined by the energy scale being probed. Applied here, 
this implies that the operator analysis is reliable provided -C 1. 

A potential fault in this estimate of the range of convergence occurs because higher di- 
mension operators may have anomalously large matrix elements. An example of this occurs 
for the dimension-six operators listed in Section [21 Consider the first dimension-six operator 
in eq.© 
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tfS(HlexpVH 2 ) 2 ] D ^^\h 2 \\ (ci~0(l)). (28) 



This should be compared to the leading quartic Higgs term coming from the dimension-five 
operators in eq.([7|) that contributes at O | 1 2 hih 2 ^j ■ One may see that the relative 
magnitude of the dimension-six to dimension-five contributions is O (^ 2 ^m ^ an @) • Thus, 

2 

strictly, the region of validity of the dimension- five operator analysis is 2fJ , M ^ an P ^ 1- 
However, as discussed in the next section, the new physics generating this dimension-six 
operator is different from that generating the dimension-five operator and so their coefficients 
should be uncorrelated. In this case the addition of higher dimension operators can reduce 
the fine tuning for some region in parameter space so that the analysis with dimension- 
five operators only will provide a useful upper bound even in regions where dimension-six 
contributions are significant. For this reason the region of validity of the dimension- five 
operators analysis is better described by the original m^/M* <C 1 and fJ,o/M* <C 1 condition. 
This keeps the corrections coming from operators with correlated coefficients small. To see this 
more explicitly consider the effect of the term in eq.(|28|) on fine-tuning. Writing v = c\ m^/M 2 
one has 

A « 2 -P M f29 ) 

p (l + 5 + 8v/g 2 )m 2 z dp { ' 

where, for example, p can be /Xq. The partial derivative is readily obtained from eq.Q. The 
dominant effect of the d=6 operator on fine tuning is the appearance of the effect of the 8vj g 2 
term in the denominator, reducing the fine tuning for the appropriate sign of 8v/g 2 . (note 
that c.f. eq. ()27p such a term is not generated by the dimension-five term at large tan /?). The 
effect of this reduction is sizeable. Taking, for example, mo/M* ~ 1/10 and c\ = 3, then 
8v/g 2 ~ 1/2 which is close to the numerical value of 5 entering the denominator. One sees 
that d = 6 operators can bring a reduction of A^2 relative to that of the MSSM including 
top/stop effects, of order (8v/g 2 )/{l + 5 + 8v / g 2 ) ps 30%. 

In the following numerical analysis we include only the effects of dimension-five operators. 
The convergence criterion found above gave mo/M* 1 and /Uo/M* C 1. In our following 
numerical analysis this bound is comfortably satisfied when we take mo/M*, /io/M* < 0.035, 
giving upper values Ci 2 < 0.07. 
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Figure 1: Left figure (a): the MSSM fine tuning A as a function of rrih; Right figure (b): the fine 
tuning in the MSSM with d = 5 operators in terms of rrih, with £i = (2 = 0.03. In both figures, the top 
pole mass considered is m t = 174 GeV for blue (dark blue) areas and = 171.2 for yellow (red) areas, 
respectively. Larger mt input (blue) shifts the plots towards higher rrih by 2-5 GeV. In both figures 
the parameters space scanned is: 1.5 < tan/3 < 10, 50 GeV < mo,mi2 < 1 TeV, 130 GeV < fio < 1 
TeV, -10 < A t < 10. 

3.4 Numerical results 

We are now in a position to determine the fine-tuning in the extended MSSM Higgs sector. 
We will present this as a function of the mass rrih of the lightest CP even Higgs. This is given 
by: 



l r 



m A + m| - yfw + £ 



+ Ci v 2 sin 2/3 



1 + 



m 2 A + m 2 z 



+ 



C2V 2 



(m 2 A — m 2 z ) cos 2 2/3 



(30) 



where 



w = [{m 2 A - m|) cos 2/3 + £] 2 + sin 2 2/3 (m 2 A + m|) 2 



m A 



mf + m| + £/2 + Ci« sin2/3- (1/2) ( 2 v 2 ; £ = 5m z sm 2 (3 



(31) 



Using the results of the previous section we compute the fine tuning for a sample of points 
in parameter space in the region with: 1.5 < tan/3 < 10, 50 GeV < mo,Tni2 < 1 TeV, 
130 GeV < no < 1 TeV, -10 < A t < 10 and 171.2 < m t < 174 GeV consistent with 
m t = 172.6 ± 1.4 GeV [29], and with the signs for £^2 chosen so as to reduce the fine tuning. 

The results are shown in Figures [T] to EJ Note that in these figures the structure apparent 
at small A and large rrih is probably a scanning artefact. We expect the under-dense wedge 
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Figure 2: Left figure (a): the fine tuning A as a function of m/j. A of MSSM is plotted in light 
blue (m t = 174 GeV) with an orange edge (m f = 171.2); A of MSSM with d=5 operators with 
Ci = C2 = 0.03 is plotted in dark blue (m t = 174 GeV) with a red edge (m t = 171.2). Right figure 
(b): similar to figure (a) but with Ci = C2 = 0.05. Non-zero or larger Q (dark blue and red areas) 
shift the plots to higher rrih to allow a reduced A for higher rrih. In both figures 1.5 < tan/3 < 10; 
50 GeV < m ,TOi 2 < 1 TeV, 130 GeV < ^0 < 1 TeV, -10 < A t < 10. 

shaped regions will be filled in with a more dense parameter sample. Similarly at very large 
A, corresponding to very precise relationships between parameters, there will be some points 
corresponding to high values of rrih that are not picked up by our finite parameter scan. 

Turning to our results, as a benchmark Figure 1(a) shows the EW scale fine tuning A of 
eq. (|13|) of the MSSM as a function of m^. One may see that A > 18 for values rrih > 114.4 
GeV, the current LEPII bound. Figure 1(b) shows A for the case of the MSSM with dimension- 
five operators added, with Ci = C2 = 0.03. The dominant effects in Figure [ljb) are mostly 
due to the effect of non-zero £1, which comes from the supersymmetric part of the higher 
dimensional operator. One may see a systematic shift of the allowed region to higher rrih 
which (for positive Q) is driven by an increase in the quartic Higgs coupling which appears 
in the denominator of the fine tuning measure (c.f. eq. (|27|) ). The overall result is that the 
minimum amount of fine-tuning A in the presence of d = 5 effective operators is small, of 
order A ~ 6, for rrih from 95 to 119 GeV. Therefore non-zero Q can accommodate larger 
rrih while keeping a A significantly smaller than in the MSSM. To illustrate the change more 
directly we superpose both plots in Figure [2] (a) . The effect is enhanced for larger operator 
coefficients, as may be seen in Figure [2] (b), where A is presented for Ci = C2 = 0.05, again 
shown relative to the MSSM case. One can see that in this case values of rrih > 114.4 GeV 
can have a low A ~ 6. Therefore A can be significantly reduced from the MSSM case, for a 
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Figure 3: As in FigEJbut with: left figure (a): Ci = 0.07, C2 = 0; right figure (b): (i = 0, C2 = 0.1 

similar m/j,. This conclusion is further supported by the plots in Figure [3] where other values 
for £j are considered. From all plots shown one sees that A < 10 is easily satisfied for values 
of the Higgs mass that can be as large as 130 GeV, depending on the exact values of Q. 

Note that, in the MSSM, A increases for low tan/3 (<S 10) and above the LEPII 
bound. However, c.f. eq. (|30p . the effect of the d = 5 operators is important for low tan/3 and 
in their presence A actually decreases for low tan /3. Thus the reduction in the fine tuning at 
very low tan /3 relative to the MSSM case is much more marked than that shown. 

The lower amount of fine tuning in the presence of effective operators is due to two 
effects. The first, already mentioned, is the presence of additional quartic Higgs couplings 
enhancing the denominator which determines the Higgs via v 2 = —m 2 /\ thus allowing for a 
smaller electroweak breaking scale. The second is the fact that higher dimensional operators 
add a tree level contribution to the Higgs mass, which reduces the need for large quantum 
contributions, and therefore the fine tuning. 

What is the scale of new physics needed for this reduction in fine tuning? Using eq. (J5J) we 
find that the scale of new physics is 

M* rj 2^o/Ci ~ (40 to 65) x no, Ci,2 = 0.05 to 0.03 (32) 

With hq between the EW scale and 1 TeV, this shows that large values of M* are allowed: 
M* » (5.2 to 8.45) TeV for fi = 130 GeV and Af* » (8 to 13) TeV for fi = 200 GeV. For 
larger fj,Q one obtains values of M* above the LHC reach. Finally, for £1 = 0.07 but with 
C2 = 0, one has M* 30 x (j,q and A < 10 for rrih ~ 130 GeV. Thus, the EW fine tuning is 
small A < 10 for 114 < < 130 GeV, for rather conservative values of £1,2- To relax these 
values one can use that an increase of £1 by 0.01 increases mh by 2 to 4 GeV for the same A. 
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3.5 The origin of "new physics" 

The presence of a higher dimension operator signals new physics and it is important to ask 
what this new physics can be. In the context of new renormalisable interactions it may come 
from the effects of new chiral superfields or from new gauge vector superfields. Consider chiral 
superfields first. One may readily obtain the d = 5 operator of eq.([T]) by integrating out a 
gauge singlet or a triplet [20j . Consider the case of a massive gauge singlet X with Lagrangian 



c x = d A ex ] x + 



d 2 e 



[iH x H 2 + X x X H X H 2 + ^M*X 2 



+ h.c. 



For M* 3> (J,, mo, one may use the eqs of motion to integrate out X, giving, to leading order 
in inverse powers of M*, 



L -x 



A 



2Af, 



d 2 e {H^f + h.c. 



(33) 



The supersymmetry breaking terms associated with this operator are obtained by replacing 
A — > X(S) giving the d = 5 operator of interest. Note that Cx has a similar form to that 
of the NMSSM. However in the NMSSM the singlet field has mass of order the electroweak 
breaking scale and cannot be integrated out whereas here we are taking the singlet mass to 
be much larger than the EW scale. 

However, the origin of the d = 5 operator cannot be uniquely ascribed to a gauge singlet 
field. Indeed it may equally well point to the existence of SU(2) triplets [SDJ [301 Ell E2] ^1,2,3 
of hypercharge ±1, 0. In this case a Lagrangian of the form 



Tle v T x + T ] 2 e v T 2 



I 



+ / d 2 9 



fiH 1 H 2 +M,T 1 T 2 + X 1 H 1 T 1 H 1 + X 2 H 2 T 2 H 2 



+h.c 



gives, to lowest order in 1/M*, eq. (f33|) with X 2 replaced by A1A2. More generally, one can 
generate the d = 5 operator through a combination of both gauge singlets and triplets. 
However note that the pure singlet X case has the advantage of not affecting the gauge 
couplings unification (at one-loop), which is not true for the SU(2) triplet. 

What about additional, massive, SU(2) doublets that couple to the MSSM Higgs sector? 
One may readily show that integrating them out does not generate, to lowest order in 1/M*, 
an operator of the type (f33|) . 
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There remains the possibility that the new physics is due to the effect of new massive vector 
gauge superfields. The simplest example is the case there is a new U(l)' gauge symmetry 
under which the Higgs sector is charged. This brings extra quartic contributions to the scalar 
potential that are expected to reduce the fine-tuning [121 ESI US]- Assuming the U(l)' is 
broken at M* one obtains the effective Lagrangian to leading order in inverse powers of M* 
given by 

q 1 Hle v H 1 +q 2 Hle v H 2 \ 2 

where g' is the U(l)' coupling and qip are the charges of the Higgses under U(l)' (qi+q% = 0). 
Note that, after including the associated supersymmetry breaking operators, this corresponds 
to the d = 6 effective operators [20] of eq.([3]) and that no d = 5 operators are generated. 

In summary, the requirement that the SUSY extension of the MSSM should not have 
significant fine tuning may indicate the presence of the d = 5 operator of eq.(pQ) which, in 
turn, suggests the presence of a massive gauge singlet and/or a SU(2) triplet. This is the 
simplest interpretation based on new renormalisable interactions but other, more complicated 
possibilities to generate the d = 5 operator may be possible. 

3.6 Further remarks on fine tuning 

Effective field theory approaches to the fine tuning of the electroweak scale were used before 
in models of low susy breaking scale scenarios [2] where both d = 5 and d = 6 operators were 
included. The model in [14J introduces supersymmetry breaking through coupling of MSSM 
states to a SM singlet field responsible for supersymmetry breaking. After integrating this 
field out, in addition to the d = 5 operator considered here, there are correlated contributions 
from the d = 6 operators. Using this, the authors find the fine tuning can be very small even 
for an arbitrarily high Higgs mass, provided the scale of supersymmetry breaking is less than 
500 GeV. 

How does this analysis relate to the one presented here? The examples given in [TJ] are 
found varying the ratio rh/M in the range 0.05 to 0.8 where rh is the supersymmetry breaking 
scale and M is the messenger mass. For rh/M small, the fine tuning is close to that in the 
MSSM but reduces rapidly for rh/M large; in this latter case the fine tuning actually reduces 
as the Higgs mass increases. This range of values for rh/M corresponds to a choice of our 
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mo/M* and [iq/M* in a similar range. The upper value strongly violates our criterion for 
applicability of the operator analysis and is a factor of ~ 10 larger than the value chosen in 
Figure [3]^ a) . Ignoring, for the moment, the fact that the contributions of higher dimension 
operators are expected to be large for this choice of mediator mass, we can ask what this 
choice of mediator mass in our analysis would give for the Higgs mass consistent with small 
A. Since the change in our upper bound on the Higgs mass roughly scales with the coefficient 
of the d = 5 operator (eq. (|30p ). this would allow a Higgs mass in the region of 276 GeV, 
much larger than our earlier conservative estimates. However, as we have stressed, for this 
value of the messenger mass the operator analysis breaks down and one should do the analysis 
including the messenger fields explicitly. 

4 Conclusions 

The LEPII lower bound on the Higgs mass places MSSM Higgs physics at the forefront of 
supersymmetry phenomenology. While this bound can be satisfied by including the MSSM 
quantum corrections, it (re)introduces some amount of fine tuning in the model. To reduce 
the fine tuning may require new physics beyond the MSSM which can be parametrised by 
higher dimensional operators. In this paper we used an effective field theory framework with 
d = 5, 6 operators in the MSSM Higgs sector, and presented a model independent approach 
to the fine tuning problem. 

We obtained exact analytical results for the EW scale fine tuning in the MSSM with 
dimension-five operators, which are also applicable to the pure MSSM case in the limit the 
coefficients of the higher dimension operators vanish. This calculation included one-loop 
corrections to the soft masses and dominant top Yukawa effects on the quartic terms of the 
potential. Similar analytical results were given for a general two-Higgs doublet model. 

Fine tuning proves to be very sensitive to the addition of higher dimensional operators and 
this is mostly due to extra corrections to the quartic couplings of the Higgs field. For the case 
of dimension-five operators we showed that one can maintain a reduced fine-tuning A < 10 
for a Higgs mass above the LEPII bound and as large as nrih ~ 130 GeV, for the parameter 
space considered, with low tan/3 (tan/3 < 10). The scale of new physics M* responsible for 
the reduction in fine tuning can be rather large, for example ~ 2/^o/Ci ~ (40 to 65) x fj,Q, 
for = 0.05 to 0.03, and M* « 30 x for (j = 0.07, (2 = 0. For values of /^o between the 
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electroweak scale and 1 TeV, these results show that large values of M* are allowed; in the 
former case M* « (5.2 to 8.45) TeV for = 130 GeV and M* » (8 to 13) TeV for = 200 
GeV. For larger /j,q, larger values of M* are possible, even above the LHC reach. These results 
follow from rather conservative choices for the coefficients of the quartic couplings induced by 
the dimension-five operators, to ensure the convergence of the effective operator expansion. 

Our numerical analysis included the effect of dimension-five operators only. These give 
the leading corrections at low tan/3, being proportional to 1/M*. However, dimension-six 
operators, suppressed by m^/M^ or fj^/M^, give contributions that can be enhanced by large 
tan/3; for (tan /3 mo) / > 1 or (tan ft fio)/M t > 1 these will be the leading terms. For this 
reason they should be included at large tan/3 and we hope to extend our analysis to the 
dimension-six case in the future. 

Of course the crucial question is what is the origin of the physics beyond the MSSM 
giving rise to these operators? The dimension-five operator can be generated by a gauge 
singlet superfield or a SU(2) triplet superfield of mass of 0(M m ) coupling to the Higgs sector. 
The dimension-six operators can be generated, for example, by an extra gauge symmetry with 
a massive gauge supermultiplet or additional (Higgs-like) SU(2) doublet supermultiplets of 
mass 0{M if ). If the fine tuning criterion is indeed of physical relevance, the significant amount 
of fine tuning found in the MSSM already indicates the need for such additional degrees of 
freedom. 
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Appendix 



A Fine tuning expressions in the MSSM. 

The coefficients Oi used in the text, Section [3.11 are: 

o-i = 0.532, a 2 (t z ) = 0.282(4.127/i t 2 - 2.783)(1.310 - hf) 1/4 

a 3 (t z ) = -0.501/i 4 2 (1.310 - /i 2 ) 1/4 , a A {t z ) = 0.532 - 5.233 h\ + 1.569/i 4 

a 5 (t z ) = 0.125 h\ (10.852 h\ - 14.221), a 6 (t z ) = -0.027 h\ (10.852 h\ - 14.221) 
a 7 (t z ) = 1-1.145/1?, a 8 (t z ) = 1.314(1.310 -h 2 t ) l/A (A-l) 

where /it is evaluated at mz and m 4 = ht(t mt ) (v / y/2) sin /3 . 



In the MSSM one obtains the following analytical expressions for fine-tuning (these are 
obtained from the results in Section [XT] by setting £i 2 = 0): 



^{(sV/8)7i sin4/3[5-(4 + 5) cos 2/3 ] + 2 [2 M g of - 7l sin2^] 74 } (A-2) 



1 



2^ 



2 (71 -/xom^c^) cos 2/3+ 73 sin 2/3 5g 2 cos /3 sin 3 /3 - (g 2 /2) sin4/3 



+ 874 



2m-o - 73 sin 2 ft + (mi2 Mo o"2 - 7l) sin 2/3 | 



(A-3) 



A, 



2 2 

mi2 r g v r 



v 2 D 



■ sin 4/3 



+2 



2^o°"2 cos 2/3 — (At (T5 mo + 2 mi2 (o"4 — oi)) sin 2/3 5 cos /3 sin 3 /3 
2m 12 o\ - ijlq02 sin 2/3 + (Af a 5 m + 2 mi 2 (cr 4 - <ti)) sin 2 /3 74 1 (A-4) 



A 



At 



|2m sin/3 [2/i o" 3 cos/3 + (2A t er 6 m - 05 mi 2 ) sin /3] (-74) (A-5) 
+ [ 5 cos /3 sin 3 /3 — - sin 4/3 ] [^o °3 c °s 2/3 — 05 m\2/2 sin 2/3 + A t mo sin 2/3] mo <? 2 v 2 j 
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and 



A £„ 2 = " 



2 B m Ho a 8 (,<■ 2 2 



v 2 D 



2 2 

\[5g 2 v 2 /l6 - 72) sin 4/3 - (4 + S)^- sin 3 2/3 - 4 7i sin 2 2/3} (A-6) 



The denominator D is now 



D = \g 2 { -5 2 w 2 (<5cos/3sin 3 /3- 1/2 sin4/3) 2 - 2 (cos 2 2/3 + <5sin 4 /3)(-7 4 )} (A-7) 



with the notation 



7i 

72 
73 
74 

and finally 

1 



/_t (-Bo m 08 + mi2 cr 2 + A t m 03) 

(-1 + (77 - A 2 CT 6 ) TOq + ^ (T 5 m 777-12 + 777, 2 2 (o"4 - <7l) +5 # 2 7J 2 /16 

2(1 - a 7 + A 2 cr 6 ) mo - -4* °"5 "ii2 m-o 

4 7 i sin 2/3 + 2 72 cos 2/3 - (4 + 5) (g 2 v 2 /8) cos 4/3 



(A-8) 



m /j,oa 8 {2 



1 



(rn 2 + 7T7 2 ) sin 2/3 



1 <^78 . 2 
1 + h> sm 2 /3 

m| + 7772 



- A«o 777120-2 - /Jo ra A t ct 3 } (A-9) 



The results for fine tuning given above considered a common bare gaugino mass, but this 
restriction can easily be lifted to obtain similar expressions. 



B Evaluation of fine-tuning A p in general two-Higgs doublet models. 

We present here the analytical result for the EW fine-tuning wrt a parameter p, for an arbitrary 
two-Higgs doublet model. This can be immediately applied to a specific model. Start with 
the general potential 



V = 7n 2 |iii| 2 + m 2 2 \H 2 \ 2 - (m§ifi- H 2 + h.c.) 



+ | A! I^Ti | 4 + | A 2 |^ 2 | 4 + A 3 | 2 



Ho 



+ A4 I Hi • H2 



+ 



- A 5 (Hi ■ H 2 ) 2 + A 6 I Hi | 2 (Hi ■ H 2 ) + X 7 \H 2 \ 2 (Hi ■ H 2 ) + h.c. 



(B-l) 



20 



In the particular case of MSSM with d = 5 operators 



^ = \{9 2 2 +9 2 l), A 2 = ±(gl+gl)(l + 5), \ 3 = ^g 2 -g 2 ) 

M = ~\gl, A 5 = C 2 , A 6 = A 7 = Ci (B-2) 

while in the MSSM alone one also sets Ci = C2 = 0. 
The minimum conditions can be written 



9 m? , dm 2 9 d\ . . 

"=-• 2X w-^w-° (B - 3) 



with 



2 ~ 2 2 , ~ 2 2 2 ~ 2 2 2 , 1^-2 ~ 2\ , /o/i / 3\ 

m = m x + m 2 — m 3 S2/3 = "^2 m 3 H 2 ' m i ~~ m i> ^\^l u ) 

A = y 4 + y 4 + (A 3 + A4 + A 5 )44 + 2A 64 S /3 + 2A 7C/34 

= v + - A 7 + 4( A 3 + A 4 + A 5 -A 2 ) + 0(l/n 3 ), (B-4) 
2 u u z 



with sp = sin/3, = cos (3, u = tan f} = vz/vx, hi = l/y/2 (v{ + hi), m 2 z = (gf + g 2 ) v 2 /4. At 
large tan (5: 

Irn 2 1 

- v 2 = — ± + — 1 (-4 m|A 2 - 8 A 7 m|) + 0(l/u 2 )- (B-5) 
A2 u A2 

Definition (|14p obtained using (|B-3[) can be used to find the most general result A p for the 
EW fine-tuning wrt a parameter p. This takes account of the dependence f3 = f3(p) induced 
by the min conditions. One finds the general expression: 

A _ dlnv 2 _ - {2w' 1 z l -(l/4)z[w 2 + [w' 3 + (l/v 2 )z> 2 } [z 3 + v 2 tgj } 
p <91nj> — (1/32) i> 2 — W3 [— z% — W4 v 2 ] 

with the following notations: 
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3 

w'i = Ag cos 4 f3 + A^4 51 cos 3 P sin p — — (A 6 — A' 7 ) sin 2 2/3 — A' 3452 cos P sin 3 /3 — A' 7 sin 4 /3, 
^2 = 4 (A 6 + A 7 ) cos 2/3+4 (A 6 - A 7 ) cos 4/3 - 2[Ai - A 2 + (Ai + A 2 - 2 A 345 ) cos 2/3] sin 2/3 
w 3 = - Ai cos 4 /3 + 2A 6 cos 3 /3 sin /3 + ^345 sin 2 2/3 + 2A 7 cos /3 sin 3 /3 + ^A 2 sin 4 /3 

w 3 = i Ai cos 4 /3 + 2A' 6 cos 3 /3sin/3 + iA 3 4 5 sin 2 2/3 + 2A 7 cos/3sin 3 /3 + iA 2 sin 4 /3 

w 4 = -(Ai + A 2 - 2A345) cos 4/3 - 2(A 6 + A 7 ) sin 2/3 + 4(A 7 - A 6 ) sin 4/3 

= — 2m 2 cos 2/3 + (m| — m 2 ) sin 2/3, 

^ = -2(m§)'cos2/3 + [(m|)' - (ml)'} sin 2/3 

4 = (m 2 )'cos 2 /3 + (m^)'sin 2 /3- (ml)' sin 2/3, 

z 3 = [4m|-4m 2 + (A 2 -Ai)f 2 ] cos2/3 + 8ml sin2/3 (B-7) 
where A 345 = A 3 + A 4 + A 5 ; A 345 j = A 3 + A 4 + A 5 - Xj, (j = 1, 2) and with: 

(»?>' - Jg, K 2 )' - Jg, - i.*>; K = ^, * = L2...7. (B-s) 

The general result ()B-6j) . (|B-7p can be applied to any two-Higgs doublet model, which includes 
all radiative corrections in the couplings and soft masses. The result in (|B-6p simplifies 
considerably in most cases, since usually many Aj are independent of p, i.e. have A^ = 0. 
It is worth taking some particular limits of the above result for A p . At large tan/3: 

1 / v 2 [ 4 A 7 (m 2 )' - 4A' 7 m 2 + A' 2 ( 2 m\ - 2 m 2 + A 3452 u 2 )] 



A p - -2 



-2 A 2 (m 2 - m 2 ,) + [-A 2 A 3452 + 2A 2 . ] t; 2 



2 [2 (mj - m 2 ) + A 3452 v 2 ] (mf/ t 
+ -2 A 2 (m 2 - m 2 ) + [-A 2 A 3452 + 2A 2 ] + tm ^ (B " 9) 



which for A 2 = A' 7 = gives 



_ 2 f _ , 2X 7 v 2 \ X 2 (m 2 )' + X 7 (m 2 )'] 

A 2 v z [ 2A 7 - A 2 [A 3452 v l + 2 (raf - m|)J 1 



In MSSM A 7 = A' 7 = then 

-2 
A 2 v 



A p = — — (m 2 2 )' + 0(1/ tan p) (B-ll) 
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which is consistent with (|B-5P and the definition of A p (assuming 6' = 0). 
In MSSM with d = 5 operators, A7 = £i so 



v 2 (1 + 5) m| 



( m 2) ^2 I ^2 /1 , x^ ( m 3i 



I - mf + m|(l + 8/2) 



+ 0(1/ tan /3) (B-12) 



which recovers the results of ([2 
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